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ABSTRACT 

Oscillating stars in binary systems are among the most interesting stellar laboratories, as these can 
provide information on the stellar parameters and stellar internal structures. Here we present a red 
giant with solar-like oscillations in an eclipsing binary observed with the NASA Kepler satellite. We 
compute stellar parameters of the red giant from spectra and the asteroseismic mass and radius from 
the oscillations. Although only one eclipse has been observed so far, we can already determine that the 
secondary is a main-sequence F star in an eccentric orbit with a semi-major axis larger than 0.5 AU 
and orbital period longer than 75 days. 

Subject headings: binaries: eclipsing - stars: oscillations - stars: individual (KIC8410637) 



1. INTRODUCTION 

The prospects for asteroseismology of red-giant stars 
have increased significan tly with the launch of Ke- 
vler dBorucki et all 120091 ) in March 2009 and CoRoT 
(jBaglin et all 120061 ) in December 2006. These satellites 
provide uninterrupted long time series of high-precision 
photometry for a large sample of stars. Among these 
stars are many red giants, allowing for both statistical 
analyses as well as detailed individual studies. CoRoT 
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data have already revealed the presence of nonradia l 
oscillation modes in red giants (|De Ridder et al.1 120091 ) 
and have given r ise to a population synthesis study 
(|Miglio et al. 2009). The currently available Kepler data 
reveal clear solar-like high radial overtone p-mode oscil- 
lations in a large sample of red giants, extending in lu- 
minosity from_th£^dchjnrp to the bottom of the giant 
branch (jBedding et all 120091 ). 

For in-depth studies, including detailed stellar model- 
ing, accurate stellar parameters are needed. Red giants 
of different masses occupy a narrow region in the H-R 
diagram, which results in relatively large errors when es- 
timating the mass from measured classical parameters, 
such as effective temperature, surface gravity and lumi- 
nosity, and from evolutionary tracks. A better estimate 
can be obtained from the characteristic frequencies of 
solar-like oscillations in these stars and their frequency 
separations, e ither through modeling o r through scal- 
ing relations dKieldsen fe Beddind IT995L Here we use 
a modeling approach and scaling relations to compute 
an asteroseismic mass and radius. 

Eclipsing binary systems provide the primary method 
for measuring masses and radii of stars. However, the 
probability of observing an eclipsing system is low, be- 
cause of the geometrical restriction that the system has 
to be viewed close to the orbital plane. We report 
here the discovery of a pulsating red giant in an eclips- 
ing binary system observed by Kepler. This star was 
found using automated variability classification methods 
to the targets of the asteroseismology progra m of Kepler 
(jDebosscher et alJ[2009l : iBlomme et afll200"9l b 

2. OBSERVATIONS 

Photometric data for KIC841063f3 (TYC 3130-2385) 
have been taken with Kepler during the commissioning 
run (Q0) from 2009 May 2 through 11, and the follow- 
ing first science run (Ql) of 33.5 d. We present results 
based on these combined time series of about 40 d. This 

19 KIC = Kepler Input Catalog. iBrown et al.l ([2005T) . 
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TABLE 1 

Stellar parameters of KIC8410637. 



KIC 



Toff [K] 
log(g) (cg.s.) 
[Fe/H] [dex] 
v sin i [km s ] 
gmicro [kms" 1 ] 



4680 ± 150 
2.8 ± 0.3 
-0.3 ± 0.3 



4650 
-0.38 



spectroscopy 
4650 ± 80 

2.70 ± 0.15 
0.0 ± 0.1 
5.0 ± 0.7 
1.3 ± 0.1 



star has a magnitude of 10.77 mag in the Kepler band- 
pass and was selected by the Kepler Asterose ismic Sci- 
ence Consortium (KASC, iGilliland et"al]|2009f ) for long- 
cadence observations. The time series data have a 29.4- 
min cadence and show both solar-like oscillations in the 
frequency regime expected for a red giant, and an eclipse 
(see Fig. [T|) . 

This star has also been part of the T rES-Lyrl ground- 
based survej0 (|Q'Donovan et al.ll2006f ) from which time 
series photometry with a time span of 75 d are at our 
disposal. In these data no eclipse is visible, while the 
primary and possibly the secondary eclipse (depth pre- 
dicted in Section 6) would have been detected, had they 
occurred during the time span of the observations. The 
fact that the primary is not seen in the TrES data implies 
that the period is longer than the duration of the TrES 
observations. 

This star was also observed photometr ically by the Su- 
perWASP cameras (|Pollacco et al.ll2006l ) for 100 days in 
2007 and for 90 days in 2008. There is a suspicion of a 
secondary eclipse, but no primary eclipse has been ob- 
served during these periods, which have some gaps due 
to weather conditions. Observations for KIC8410637 
are also ava ilable from the All Sky Automated Sur- 
vey (AS AS) (jPoima nski 1993), in which there are three 
points that have lower flux than the majority. These 
points might be due to the eclipse, which would imply 
an orbital period of about 380 d or a submultiple thereof. 
Unfortunately, there are too few points, and their uncer- 
tainty is too large to draw firm conclusions. 

Additionally, we checked whether we could put more 
constraints on the period from the time stamps of the 
data we have at our disposal. We determined which peri- 
ods would have led us to miss primary eclipses because of 
a lack of observations. We would have missed all eclipses 
for an orbital period of about 135, 220, 253, 261, 270 and 
295 days. All other possible periods are longer than the 
380 days, mentioned above. 

3. STELLAR PARAMETERS 

For KIC8 4 10637 the Kepler Input Catalogue 
(jBrown et al.l 120051 ) provides a set of s tellar pa- 
rame ters, as listed in Table ffl Tycho- 2 (|H0g et al.l 
120001 ) and 2MASS dSkrutskie et all (20061 ) photometry 
are also available. lOfekl (|2008l ) fitted the se colors wit h 
synthetic photometr y of stellar t emplates ( Pickles 1998T ). 
Using two methods. lOfekl (|2008l ) classified the star as a 
metal-weak K2III star with an effective temperature of 
4650 K or a K4V star with a temperature of 4350 K. 

We also ob tained three spectra with th e HERMES 
spectrograph (jRaskin fe Van Wh ickcl 20Ql) mounted on 
the Mercator Telescope, La Palma, Spain: two on 2009 

20 http: //nsted. ipac . caltech. edu/NStED/docs/holdings .html 
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Fig. 1. — Top: time series data from Kepler of KIC8410637. 
Middle: a section of the light curve centered on the eclipse with the 
best-fit model (solid line). Bottom: the residuals after correcting 
for the eclipse model, clearly showing the solar-like oscillations. 

October 13 and one on 2009 October 15. These spec- 
tra have a resolution of about 85 000 and a wavelength 
rang e of 4000 - 9000 A. We analyzed the average spec- 
trum, which has a typical signal-to-noise ratio of ~ 90 at 
6000 A , with methods described by iHekker fc Melendezl 
(|2007l ) and list the results in the last column of Ta- 
ble [T] This method does not directly provide uncer- 
tainties on the individual parameters. Here, we adopted 
the tot al of the difference a n d sca tter between the val- 
ues of IHekker fc Melendezl (120071) and the values of 
ILuck fc Heiteii (|2007l ) (see IHekker fc MelendedllM for 
more details). For wsini, we computed the uncertainty 
from a 10% uncertainty in the total FWHM of the lines, 
which is combined with the a on the macro turbulence. 
These results are in agreement with those in the KIC and 
the literature. They confirm the red-giant nature of the 
primary star. 

We investigated the spectrum for signatures of the sec- 
ondary star, in the form of additional absorption lines 
and cross-correlation profiles, but no clear signatures 
could be detected. The three spectra obtained over a 
time span of 3 days are not sufficient to see radial veloc- 
ity variations due to the binary orbit. 

4. ASTEROSEISMIC ANALYSIS 

Solar-like oscillations are visible both during and out- 
side the eclipse (Fig. Q]). To analyze these we adopted 
two approaches. In the first approach, we corrected the 
light curve including the eclipse by fitting and subtract- 
ing trends for different parts of the time series using a 
linear polynomial for the commissioning data, a second- 
order polynomial for the part of the Ql data before the 
eclipse, the mean value for the data in full eclipse, and a 
linear polynomial for the data obtained after the eclipse. 
We omitted the observations during ingress and egress. 
The resulting time series is shown in the top panel of 
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Fig. 2. — Top: flux as a function of time, after removing the 
binary signature as described in the text. Centre: the power spec- 
trum of the red giant oscillations. The power excess due to solar- 
like oscillations is present in the range 30 — 60 fiHz. The excess 
below 20 /iHz might be due to granulation in the primary's at- 
mosp here. Bottom: H-R diagram with BASTI iPictrinfcr ni et aLl 
I2004T) evolutionary tracks for different masses and solar metallicity. 
The filled circle shows the position of the red giant. The shaded 
area shows the secondary position constrained by two values de- 
rived from the light-curve fit: the radius ratio (dashed lines) and 
the bolometric luminosity ratio (dashed-dotted lines; taking into 
account the Kepler bandpass). The asterisk indicates the position 
of the Sun. 

Fig. [51 In the second approach, we omitted the complete 
eclipse from the time series. 

The Fourier spectrum shows a clear power excess (see 
the second panel of Fig. [5]). The power excess is present 
between 30 — 60 /iHz, which is the dominant range of 
^max (frequency of maximum oscillation power) found 
for a large sample of red giants observed by CoRoT 
(jHekker et al.ll2009bl ). This reflects the high population 
density of red-cl ump stars, compar ed to stars on the as- 
cending branch (|Miglio et al.ll2009f ). No power excess is 
present at frequencies between 70 /xHz and the Nyquist 
frequency of ~ 280 fiRz. We analyzed the power spec- 
trum obtained from the correc ted time series with two 
methods of the Octave pipeline (jHekker et al.ll20"09al : Oc- 
tavel and Octave2) and the power spectrum of the light 
curve where the eclipse had been om itted with the auto - 
mated pipeline developed in Sydney (|Huber et al.ll2009h . 

Octavcl computed ^ max and A^—q (maximum mode 
amplitude) from a binned power spectrum, and Av (large 
frequency separation) using the power spectrum of the 



power spectrum. Octave2 computed v max and Ai = q 
from a Gaussian fitting to the oscillation envelope, and 
Av from an autocorrelation of the oscillation frequen- 
cies. The latter have been determined with a Baycsian 
method, where oscillation frequencies are those with a 
posterior probability of less than 0.1% of being noise. 
The Sydney pipeline produced ^ max from a smoothed 
power spectrum. The current data do not yet allow us to 
measure the life times of the solar-like osc illations. For 
more d etails on the pipe l ines w e refer to iHekker et alJ 
(|2009ah and lHuber et all (f2009h . 

The resulting parameters are listed in Table [H Note 
that Ai = q is computed for a power spectrum with mean 
squared scaling, which is a facto r y/2 different with 
respec t to the scaling relations of iKieldsen fc Beddinel 
(| 19951 ) . With the effective temperatu re and the scaling 
laws from Kjeldsen & Bedding (1995) wc computed the 
radius, mass and luminosity of the primary (see Table[5]). 

We also used the RADIUS pipeline described by 
iStello et al.l (|2009f ) to compute the stellar radius. The in- 
put parameters for this pipeline are the large separation, 
T e ff and metallicity, for which we used the spectroscopic 
values listed in Table [T] The results from this analysis 
are listed in Tableland lead to a value for log(g) of 2.6 
±0.1, which is in agreement with our spectroscopically 
derived value. 

The values computed with different methods and 
approaches are in good agreement. The weighted 
mean values of the results from the different meth- 
ods are used throughout the rest of the paper, 
i.e., R R = 11.8 ± 0.6 R Q , L K = 58 ± 6 L , 
Mr = 1.7 ± 0.3 M Q (suffix 'R' refers to the red gi- 
ant) together with the spectroscopically determined ef- 
fective temperature. We point out that the estimates of 
these fundamental parameters rely on the input physics 
of stellar models, through the interpretation of obscrv- 
ables (e.g. v max and Ai = q) in terms of the properties of 
the convection and thermodynamics, which we assumed 
to be known and error-free. This implies that there may 
be unknown systematic uncertainties. Nevertheless, the 
radius and mass have values in good agreement with the 
expected range for red-clump stars. For the radius there 
is also agreement between the result from scaling rela- 
tions and from models (RADIUS). 

5. BINARY ANALYSIS 

The very flat bottom of the eclipse, except for the red- 
giant oscillations, seems to be indicative of a smaller ob- 
ject occulted by a considerably larger object. This is 
underlined by the fact that for an annular eclipse, i.e., 
transit, limb darkening effects play a more prominent role 
and would have caused a more gradually changing shape 
of the eclipse. Because wc have only one eclipse, a full 
analysis of the binary is not yet possible, but we are able 
to put several constraints on the nature of the secondary 
and on the orbit. 

From the Kepler data, we directly measured the depth 
of the eclipse to be about 9%. This depth is determined 
by the luminosity ratio of the two stars, i.e., I = LijL^ 
= 0.099 ± 0.005. The uncertainty is estimated taking the 
shape of the light curve outside the eclipse into account. 
In addition to the depth, wc measured the time dur- 
ing which the secondary is fully obscured: 1.64 ± 0.01 d. 
The ingress and egress are symmetric and each lasts for 
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TABLE 2 

Oscillation parameters of KIC8410637 computed with different methods. 





Octavcl 


Octave2 
measured 


Sydney 


RADIUS 
model 




45.0 ± 4.8 


45.2 ± 1.3 


44.2 ± 0.5 




Ai/ [/iHz] 


4.64 ± 0.23 


4.5 ± 0.1 


4.51 ± 0.37 




8i/ 02 [/LtHz] 


0.6 ± 0.1 




0.69 ± 0.09 




Af =0 [ppm] 


57.6 ± 7.3 


61.8 ± 6.2 


53.4 ± 1.9 








derived 






fl [Rq] 


11.2 ± 1.6 


11.8 ± 0.7 


11.4 ± 1.9 


13.1 ± 1.7 


L [L ] 


53 ± 16 


58 ± 8 


55 ± 19 


72 ± 19 


M [Mq] 


1.7 ± 0.7 


1.8 ± 0.3 


1.7 ± 0.9 





0. 276 ± 0.007 d. The symmetry of the ingress and egress 
indicates that the relative acceleration of the binary com- 
ponents between ingress and egress is equal to zero, i.e. it 
is justified to assume the sky-projected relative velocity 
of the two stars to be constant during the eclipse. 

The eclipse time depends on the radii of both stars and 
the geometry of the system. Knowing that the eclipse 
lasts only during ~ 2% or less of the total orbital period, 

1. e., only a very small segment of the orbit is covered 
during the eclipse, and that the system is viewed close 
to the orbital plane, implies that the time of total eclipse 
(Ttotai) can be expressed as: 



Ttotal 



P(R R cosS- R 2 ) 
AaE{e) 



(1) 



with P the orbital period, a the semi-major axis, E the 
complete elliptic integral, e the eccentricity and 5 the 
latitude of the eclipse on the stellar disc and R R and R 2 
the radius of the red giant and the secondary, respec- 
tively. The total time of ingress, total eclipse and egress 
(■Hn+totai+eg) from the same geometric considerations can 
be expressed as: 



'Hn+total+cg 



P(R R cos S + R 2 ) 



(2) 



4aE(e) 

Dividing Eq. (fT]) by Eq. provides an estimate of R 2 
as a function of cos S: 

T to t a i Rr cos S - R 2 cos 5 - r 

= x = = , (3) 

Tin+total+eg Rr COS 5 + R 2 COS + T 

with r = R 2 / Rr,. In this way we found that 

R 2 ,1-x 

— = cos , 

R R 1 + x' 



(4) 



where x = 0.75 ± 0.03. Using 6 = (central eclipse) 
and the value of the radius of the red-giant star derived 
in the previous section, we then find an upper limit of 
R 2 < 1.7±0.1R Q . 

To confirm this geometric analysis we computed a ba- 
sic model of the eclipse signal. We neglected limb dark- 
ening and computed the fractional overlap between the 
two discs, taking into account the relative luminosity of 
the star being eclipsed. As explained above, we assumed 
the sky-projected relative velocity (v) of the two stars 
to be constant. This model has been fitted to the por- 
tion of the light curve near the eclips e using a Markov- 
Chain Monte Carlo approach (see e.g. iTorres et al.ll2008l . 
and references therein), which yields the joint posterior 
probability density of all the parameters, given the data. 



The physical parameters in our model are r, v, I, and 
b = sin S (the impact parameter). In order to avoid 
strong correlations among these parameters (in partic- 
ular between r, v, and b), when performing the fit we 
replace v with u 2 /(l — b 2 ) and r with r/v 2 . The values 
of the physical parameters resulting from this fit are b = 



0.0+° 



.15 



0.1435 



+0.0008 
0.0027' 



1.0338 



+0.0005 
-0.011 



-%/day, 



and I = 0.10060lo;°ooo6> consistent with the values ob- 
tained from the geometric analysis. 

To test our assumption that v is constant, we repeated 
the fit with allowing v to change linearly with time. The 
result was a deceleration of 0.015t^^i? R /day 2 (i.e. ~ 
3% slower at egress than at ingress) . This may be a real 
asymmetry in the eclipse, implying an eccentric orbit. 
However, it could also be due to the pulsation signal 
altering the apparent shapes of the ingress and egress. 

We converted the fitted value of I into a bolometric 
luminosity ratio by folding blackbody curves with the 
Kepler spectral response curve. This is a function of the 
secondary's temperature, as shown in Fig. [5] (dot-dashed 
lines). With the constraint from the radius ratio and 



£2 
Lr 



ft 2_ t effj 
R2 7.4 > 
-"-R /eff,R 



(5) 



we find L 2 = 5.2 ± 0.7 L Q . The secondary is thus proba- 
bly an F main-sequence star. Using the mass-luminosity 
relation for stars with masses between 0.5 and 2.0 M Q , 
i.e., (L/Lq) = (M/Mq) 4 - 5 , we found a mass for the sec- 
ondary M 2 of 1.44 ± 0.05 Mq, where we did not take the 
uncertainty of the mass- luminosity relation into account. 
We also computed the surface gravity for the secondary 
and found log(g 2 ) = 4.2 ± 0.1. 
The binary also has to obey Kepler's third law: 



M R + M 2 



4ir 2 a 



2„3 



(6) 



GP 2 ' 

in which we still have a as unknown, and we take a min- 
imum value of 75 d for P. With this value we computed 
a minimum semi-major axis for the system of 0.5 AU. 

On the other hand, if we assume a circular orbit, we 
can compute the period. From the eclipse fit we know the 
orbital velocity and combining this with Kepler's third 
law, we get: 

2ttG(M r + M 2 ) (7) 



P 



Using the masses of both stars, we obtain P = 32 ± 6 
days. This is considerably shorter than the minimum 
orbital period we inferred from currently available obser- 
vations, which implies an eccentric orbit. 
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TABLE 3 

PARAMETER OVERVIEW OF ECLIPSING BINARY KIC8410637. 





primary 


secondary 


R [R Q ] 


11.8 ± 0.6 


1.7 ± 0.1 


L [L ] 


58 ± 6 


5.2 ± 0.7 


M [M Q ] 


1.7 ± 0.3 


1.44 ± 0.05 


T eS [K] 


4650 ± 80 


6700 ± 200 


log(g) (cg.s.) 


2.70 ± 0.15 


4.2 ± 0.1 


orbital parameters 


-Porbit [days] 


> 75 




a [AU] 


> 0.5 





6. DISCUSSION 

From observations obtained with Kepler, we found a 
pulsating red giant in an eclipsing binary. Once a full 
orbit has been observed, the orbital parameters will pro- 
vide an independent measure of the stellar parameters, 
such as mass and radius, with respect to the asteroseis- 
mic values. These additional constraints make this a very 
interesting case and the star is therefore being followed 
up by the Kepler Mission as well as with ground-based 
spectroscopy. 

Stellar parameters of the primary red giant have been 
computed from a spectroscopic analysis and from the 
solar-like oscillations. We did not find evidence of the 
secondary component in the spectra. The currently avail- 
able spectra are not yet sufficient to put constraints on 
the orbit. 

From the single eclipse observed by Kepler and addi- 
tional independent observations from TrES, SuperWASP 
and AS AS, wc inferred constraints on the orbit and sec- 
ondary star. The annular eclipse, i.e., the secondary 
passing in front of the red giant, would cause an eclipse 
with a depth of ~ 1.7%, assuming a circular orbit and 
without taking limb-darkening into account. Such a dip 
would be observable, but has not been detected by Ke- 
pler, TrES or AS AS. Only the SuperWASP data show 
a feature that might be due to an annular eclipse. The 
SuperWASP and ASAS data have gaps due to weather 
conditions. Although there are also three gaps (4.7, 2.9 
and 1.7 d) in the TrES time series data with a duration 
longer than the duration of the eclipse, we expect the 
chance that an eclipse falls exactly during one of these 
gaps to be low. Therefore, we infer that the orbit is 
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longer than the time span of the TrES data, i.e., at least 
75 d. Also, an initial inspection of the raw Kepler data of 
Q2 (additional 3 months of data following the Ql phase) 
does not seem to reveal another eclipse. In the case that 
we missed an eclipse due to a lack of observations, the 
minimum orbital period would be 135 days. 

From the orbital velocity of the secondary, we com- 
puted the period for a circular orbit, which would be 
shorter than half the total time span of the TrES data 
and we should have seen an occultation in those data. 
Therefore we concluded that the orbit is eccentric. The 
possible annular eclipse in the SuperWASP data has a 
length of ~ 7 days and depth of ~ 0.03 mag. If confirmed 
this would indeed imply a large orbital eccentricity. 

From the eclipse times and depth we were able to com- 
pute the luminosity and the radius of the secondary. All 
values seem to be compatible with an F main-sequence 
star. The mass-luminosity relation for main-sequence 
stars then leads to a mass estimate of the secondary as 
well as its surface gravity. These results are summarized 
in Tableland the position of both stars are indicated in 
an H-R diagram in Fig. [5J 

KIC8410637 is an extremely interesting binary for fur- 
ther follow-up. Longer time series from Kepler will im- 
prove dramatically the detection threshold for the deriva- 
tion of the properties of the oscillations of the red giant. 
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